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Simple patterns for non-linear susceptibilities near Tc 
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Non-linear susceptibilities upto the eighth order have been constructed in QCD with 2 flavours of 
dynamical quarks. Beyond leading order, they exhibit peaks at the cross over temperature, Tc. By 
analyzing their behaviour in detail, we find that the dominant contributions near Tc come from a 
set of operators with a remarkably simple topology. Any effective theory of QCD near Tc must be 
able to explain these regularities. 
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Quark number susceptibilities (QNS) in QCD P| are interesting because they are measurable through event-to-event 
fluctuations of conserved quantities in heavy-ion collisions 0- Recent determinations of the linear QNS in lattice 
QCD include those in the continuum limit of the quenched theory the first results in the high temperature phase 
of Nf = 2 QCD and the first computation in TV/ = 2 + 1 QCD [||. The non-hnear susceptibilities (NLS) are 

a generalization introduced in and have been used in finding the Taylor expansion of the pressure of the QCD 

plasma at finite chemical potential. The linear combinations used for pressure were also reported in Nf — 2 QCD 

Here we report on systematic simplicities of these quantities that we discovered in our investigation of QCD with 
light dynamical quarks. These simple patterns which we find here for the first time may be consistent with weak 
coupling theory in the high temperature phase of QCD. However, in the vicinity of the finite temperature cross over 
at Tc, we find a different simple pattern. It seems possible to incorporate it into a simple model of the physics of the 
cross over. A few of these results have been discussed in Here we complete the study of the NLS started there. 

The partition function for QCD at temperature T and chemical potentials /i/ for each of Nf flavours, can be written 
in the form 
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where Sg is the gluon part of the action and M denotes the Dirac operator. The pressure, 
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FIG. 1: X2o/T^ varies smoothly across Tc, and behaves roughly as an order parameter, being small in the hadronic phase and 
large in the plasma. X^i/T^ is small in the hadronic phase, perhaps peaks near Tc and is not significant in the plasma phase. 
Data from lattice sizes 4 x 16'^ (circles) and 4 x 24^ (boxes) are shown. 
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FIG. 2: The operator 0ii is shown on the left with the smallest number of gluon connections between the two fermion loops. 
The contribution is naively of order g®, but a computation shows that it is actually g^logg [l^ . The operator 02 shown on 
the right is fermion line connected and hence the leading contribution shown is of order 1. The black dots denote insertions of 
7o arising from the derivatives with respect to the chemical potential. 



which is a convex function of T and fij, can be expanded in a Taylor series about the point where all the /i/ = j^]. 

In this paper we examine staggered fermions with Nf = 2 and a small but non- vanishing quark mass, = md = m. 
The iV-th order derivatives in the Taylor expansion then can be taken n„ times with respect to fiu and Ud = N — 
times with respect to /id- This is the non-linear quark number susceptibility (NLS), which we write as Xnu.nd- This 
new notation streamlines a more cumbersome notation which was used earlier. The translation table between these 
two notations can be understood from the relations — 

X20 = Xuu = Xdd = X02, Xll = Xud, Xia = Xuuuu — Xdddd = X04, X22 = Xtmdd, etc, (3) 

where we have used flavour symmetry to write Xn„,nd = Xnd.n^- 
The Taylor expansion of the pressure can be written as 

AP{T, ^l^,^ld) = P{T, ^l^,^ld) - p{t, o, o) = V x«„,„, ^ (4) 

The NLS above can be written down in terms of the derivatives of Z. From the expression in eq. it is clear that 
the derivatives with respect to the fif land entirely on the determinants. Now, since Det A/ = expTr logM, the first 
derivative gives (Det M)' — Tr (A/~^Af')Det M = 0iDet M. Higher derivatives can be found systematically using the 
additional relation MM^^ = 1, which yields (M"^)' = -M"iM'M"^ Our notation for operators is that 0^ = 0„+i, 
and 0imn- - — &i&m&n ' ' •• The expectation values (02n+i(M/ = 0)) = by CP symmetry. The derivatives of Z can 
be written in terms of expectation values of certain operators involving powers of traces of products of inverses and 
derivatives of the Dirac operator. Diagrammatic methods for their evaluation were developed in 0, |^ and explicit 
expressions were written down in p^ . 

We report on results obtained using the configurations generated in the study reported in Details of our 

simulations and statistics can be found there. These results have been obtained on lattices with temporal extent 
Nt = 4, and varying iV^, with the spatial volume being large. The quark mass has been fixed in physical units to be 
such that m-^/mp — 0.31 ± 0.01, about 50% larger than in the real world, making this the smallest quark mass at 
which NLS have been studied. Details of how the temperature scale is set on the lattice can also be found in In 
physical units we find that the cross over temperature, Tc is m„/ Tc = 5.4 ± 0.2. 

The volume dependence of Tc has been remarked upon in [lO| ; we see evidence of some volume dependence in the 
bare coupling at the cross over, but the scale has larger uncertainties, so a finite size scaling study of the shift of Tc 
with V performed at these lattice cutoffs a will not be very useful. However, strong finite volume effects on the NLS 
were found when the spatial lattice extent was too small, Ng < ANt. In the remainder of this study, therefore, we 
concetrate on the NLS obtained with Ng — 16, using data obtained with Ns = 24 to make cross checks of the results. 
At Tc, the finite volume shift in the results is significant, but become negligible on moving slightly away — to 0.95rc 
or 1.05Tc, for example. 

The two leading terms in the series, the diagonal QNS, X20, and the off-diagonal QNS, xiii have been computed 
before. For completeness we display results from lOj in Figure^ Note that xii = {T/V){0ii), which is a quark-line 
disconnected diagram. Also, one can see that X20 ~ Xii = {T/V)02, which is quark-line connected. Diagrammatic 
representations of these are shown in Figure [3 Recall that for T > T^ these diagrams have been computed in weak 
coupling theory, giving reasonable agreement with the lattice results j^, 13, 14j. 

A counting rule for the minimum number of gluon lines needed in a quark-line disconnected diagram was obtained 
in [l^ by noting that effectively the diagrams are Abelian, and Furry's theorem holds, i.e., the number of 7^ insertions 
must be even. Among these must be counted the insertions of 70 arising from taking derivatives with respect to the 
chemical potential. For 0ii one gluon exchange is ruled out for reasons of gauge invariance, two by the counting rule, 
and hence three gluons are needed, as shown in Figure |21 
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FIG. 3: X40 (boxes) peaks at T^, and the peak is entirely due to the term in {T/V){022) (circles), as shown on the left. After 
subtracting this out, one gets a much smoother function (circles in the right panel), which agrees well with {T/V){04) (boxes). 
Data are from lattice sizes 4 x 16'^. 




FIG. 4: One of the contributions to the operator 022 is shown on the left with the smallest number of gluon connections 
between the two fermion loops allowed by the counting rules of [l^ . Other contributions correspond to permuting the gluon 
lines and operator insertions along each quark line, while keeping the loop topology fixed. The operator 04 shown on the right 
is fermion line connected and hence of order 1. These diagrams are expected to give accurate results in the plasma phase. 



In figure n some volume dependence is visible in the immediate vicinity of Tc- The high temperature behaviour of 
X2o/T'^ is consistent with our earlier results in 6], and, therefore, is compatible with the predictions of [T^ [T3l|. The 
results on xn/T^ are also completely compatible with earlier results in |a| after correcting for a division by an extra 
factor of (T/V) for Xii/^^^ reported there. Comparison with the recent results of 0,0 are harder to perform since 
the actions and quark masses are different. 

At the fourth order, there are five operators — 04, O31, 022, 0ii2 and 0iiii. The last four are quark- line discon- 
nected operators. The connected parts of the operators enter into the expressions for the NLS |3ill3- ^^i^ paper, 
we decompose the NLS into connected parts of these operators, such as (r/V^)(022)c- Since comparisons are always 
with connected parts, wc indulge in slight notation-abuse by dropping the subscript often. We remind the reader of 
the definitions of the connected parts at the fourth order — 




^031^ and ^^4^ connected pieces by themselves; the former by virtue of the fact that (0„) — for odd n, the 

latter because it is the largest loop at this order. In ^| we have shown that each distinct operator topology is a 
physical observable in a version of QCD with appropriate number of quark flavours. 

We show our results for the QNS in Figure 01 where we also plot the connected part of (T/V^)(022) multiplied 




FIG. 5: xii (boxes) and (T/V){@33) (circles) obtained on a 4 x 16^ lattice. 



by the coefficient witli whicli it enters into X40- Tfiis operator appears to talce care of the peak in the QNS near T^. 
In Figure|21we have also shown the difference between these two quantities. The peak disappears and the remainder, 
in the high temperature phase, is saturated by {T /V){0i) . Like 02, this expectation value is also like an order 
parameter, being small in the low T phase, and large on the other side of T^. 

In this range of temperatures, the two major contributions to the fourth order QNS are from 04 and 022- We 
also find this kind of peak in X22, where it again matched the peak in {T /V){022) ■ No other QNS at this order has 
contribution from this operator, and also show little sign of a comparable peak near Tc- 04 gives no contribution to 
any other QNS, and, compatible with this, we see that all other 4th order QNS are very small above Tc- A similar 
behaviour is also seen on the 4 x 24^ lattice. 

It is interesting that the counting rules of show that the two largest contributions above Tc should come from 
precisely these operators. 04 is of order 1, and the connected part of 022 shown in Figure 01 is naively of order g**. 
In comparison, 03i is of order 5^, 0ii2 is of order and 0iiii is of order g^^ . These naive powers may be modified 
into some logarithms in the computation. 

At the sixth order we have eleven topologically distinct operators 06, 05i, 042, 033, 0ii4, 0i23, 0222, 0iii3, 
01122, 011112 and 0111111- The determination of the NLS are also significantly more expensive than the linear QNS, 
requiring many more vectors in the stoachstic evaluation of the traces [T^ • One result is that the measurements are 
more noisy at higher orders. Nevertheless, it is possible to make significant statements about the structure of these 
operators. 

One interesting point, illustrated in Figure [S] is the qualitative similarity between xii a-nd (r/l^)(033). Both are 
small in the high T phase, possibly peak in the vicinity of Tc, and are comparable to other operators in the low T 
phase. We have previously argued that the increase in the ratio X11/X20 with decreasing T implies that the fermion 
sign problem becomes more severe, thus restricting the usefulness of all the recent methods which have been developed 
to handle this problem. The observation in Figure [51 extends this argument to finite chemical potential. 

However, at higher temperatures, such contributions are small. The dependence of xeo on T is shown in Figure 
The peak at Tc is due to contributions from 222, as we demonstrate by plotting along with this the values of X42 
normalized so that the two have equal contribution from 0222- The difference is small; for T > Tc it is saturated by 
06, which is much smaller than the peak, but much larger than 222- The operator 024 also peaks at Tc, but the 
value at the peak is negligible in comparison with 0222- Power counting shows that Og is of order 1, 024 is of order 
5'*, but 0222 is of order g^ . The form of the operators is shown in Figure [T] This is the lowest order at which we first 
find explicitly that the perturbative power counting of the high temperature phase does not extend down to Tc. 
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FIG. 6: In the first panel we show xeo (circles) and 5x42 (boxes) as found on a 4 x 16"^ lattice. The two are normalized such 
that they have equal contribution from 0222- The second panel shows {T/V){0f})T^ on 4 x 16^ (circles) and 4 x 24^ (boxes) 
lattices. Note the difference in the scales of the two figures. 




FIG. 7: One of the contributions to the operator 0222 is shown on the left with the smallest number of gluon connections 
between the two fermion loops allowed by the counting rules of fl^ . Other contributions correspond to permuting the gluon 
lines and operator insertions along each quark line, while keeping the loop topology fixed. The operator 06 shown on the right 
is fermion line connected and hence of order 1. These diagrams are expected to give accurate results in the plasma phase. 
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FIG. 8: In the first panel we show xso (circles) and 7x62 (boxes) as found on a 4 x 16'^ lattice. The two are normalized such 
that they have equal contribution from 02222. The second panel shows (r/V^)(08)r* on a 4 x 16^ lattice. Note the difference 
in the scales of the two figures. 
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FIG. 9: The connected parts of the expectation values of 026 (circles) and 044 (boxes) as found on a 4 x 16"^ lattice. The 
expectation values are normalized by T/V and rendered dimensionless through a multiplication by T*. 



This pattern recurs at the eighth order, as we display in Figure |H1 There is a peak in some of the susceptibilities 
at Tc, but this can be ascribed to 2222- The high temperature phase is dominated by a non-vanishing value of 08, 
which is much lower than the peak. Other operators at the eighth order which may peak at Tc are 026 and 044. As 
we illustrate in Figure IHl they indeed have interesting behaviour near Tc. However, these operators are numerically 
negligible compared to the value of 2222. In the high temperature phase the power counting rules show that 0s is 
of order 1, 026 and 044 are of order g'^, whereas 2222 is of order g^. The pattern of dominance near Tc therefore has 
nothing to do with power counting in g. 

In summary then, we have found a very pleasing pattern for the NLS. In the hadronic phase, all operators seem 
to have comparable expectation values. This is not unexpected. In the hadronic vacuum, at T = 0, many different 
operators have vacuum expectation values, which are all typically expected to be of similar order. Above Tc, we have 
an extremely simple pattern, in which the NLS are dominated by the operators with a single quark loop, 0„, and the 
expectation values (r/F)(0„)r"~^ are all in the range of 1-2. This pattern seems to be organized by weak-coupling 
power counting arguments, but it would be useful to have precise estimates of these operators through perturbative 
computations. 

It follows from this observation, that the pressure at finite chemical potential has contributions from all even terms, 
but the numerical importance of the terms decreases factorially at high temperature. As shown in in a free 
field theory at finite fj,, the pressure can be separated into a quark piece and an antiquark piece, each of which has 
contributions to all even orders in /x, which cancel to give a pressure which contains only terms upto order /x^. These 
small terms in the pressure can be thought of as little shift in these pieces caused by a weak coupling, such that the 
cancellation becomes incomplete. Such a mismatch between particle and antiparticle is possible because a chemical 
potential explicitly breaks CP invariance. 

The most unexpected regularity that we have found is in the vicinity of Tc. Here, the NLS are dominated by a 
composite operator which is made up of appropriate numbers of fermion loops with two 70 insertions in each, i.e., with 
an appropriate number of 02. Our observation suggests that it may be possible to write down effective long-distance 
theories in which this composite bosonic operator is treated as a field operator whose expectation value shows the 
correct cross over behaviour. In that case (022)c would be the susceptibility of this field, and being proportional to 
the temperature derivative of (02), would peak, as observed. The expectation value (0222)c would be proportional 
to the next derivative of (02). Then the T-dependence of these quantities at /i/ = would have the shapes shown 
Figure Cni 

In an effective 3-d spatial Landau theory of the form that we suggest, 02 can be taken to be a two point function built 
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FIG. 10: The expectations for the NLS near Tc based on an effective theory of QCD near the phase transition in which the 
composite operator 02 is identified as the order parameter. 



from one polarization of a vector operator. Under the symmetries of the transfer matrix that builds the equilibrium 
correlation functions, i.e., the screening correlators, this polarization mixes with the scalar 1^. It has been suggested 
that the scalar crucially impacts the physics of the phase transition in the chiral limit |l7|, because of the fact that 
it becomes massless at that point. This is the situation in the chiral limit; it would be interesting to see predictions 
from such models for the behaviour of these NLS at finite quark mass. 
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